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QED-model for multichain coupled qubit system, proposed in was confirmed by Raman scat- 
tering studies of carbon zigzag-shaped nanotubes, produced by high energy ion beam modification of 
natural diamond single crystals. New quantum optics phenomenon - Rabi waves - has been experi- 
mentally identified for the first time. Raman spectra in perfect quasi-lD carbon nanotubes are quite 
different in comparison with well known Raman spectra in 2D carbon nanotubes of larger diameter. 
They characterized by vibronic mode of Su-Schriffer-Heeger cr-polaron lattice and its revival part in 
^ . frequency representation, which is the consequence of Rabi wave packet formation. 
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I. INTRODUCTION AND BACKGROUND 



Quantum electrodynamics (QED) takes on more and more significance for its practical application and it, in 
fact, becomes to be working instrument in spectroscopy studies and industrial spectroscopy control. QED-model for 
multichain coupled qubit system was proposed in [lj. It is generalization of the model described in generalizing, 
in its turn, Tavis-Cummings model by taking into account the ID-coupling between qubits. The most substantial 
result in [2j is the prediction of new quantum optics phenomenon - Rabi waves formation. It is substantial, that 
in the model, proposed in [1| the interaction of quantized EM-field with multichain qubit system is considered by 
taking into account both the intrachain and interchain qubit coupling. Tne model is presented with an example of 
perfect zigzag-shaped carbon nanotubes (ZSCNTs). ZSCNTs can be considered to be the set of n carbon backbones 
of trans-polyacetylene (t-PA) chains, which are connected between themselves. Given set can be considered to be a 
single whole, which holds the quasionedimensionality of a single chain. It seems to be correct for perfect ZSCNTs, if 
k*" ■ their diameter is < lnm (see Sec. II). 

It is substantial, that the physical properties of t-PA chains are studied very well both theoretically and experimen- 
tally. Let us touch on given subject in more detail. The authors of Q, [H have found the most simple way to describe 
mathematically the chain of t-PA by considering it to be Fermi liquid. It is now well known SSH-model. The most 
substantial suggestion in SSH-model is concerned the number of degrees of freedom. Su, Schrieffer, Heeger suggested, 
that the only dimerization coordinate u n of the n-th C-ff-group, n — 1, N along chain molecular-symmetry axis x 
is essential for determination of main physical properties of t-PA. Other five degrees of freedom were not taken into 
consideration. Nevertheless, the model has obtained magnificent experimental confirmation. In [l[ was suggested, 
that given success is the consequence of some general principle and it was shown, that given general principle is really 
exists and main idea was proposed by Slater at the earliest stage of quantum physics era already in 1924. It is - "Any 
atom may in fact be supposed to communicate with other atoms all the time it is in stationary state, by means of 
virtual field of radiation, originating from oscillators having the frequencies of possible quantum transitions..." 
The development of given idea is based on the results of the work Q. It has been found in that Coulomb field in 
ID-systems or 2D-systems can be quantized, that is, it has the character of radiation field and it can exist without 
the sources, which have created given field. Consequently, Coulomb field can be considered to be "virtual" field in 
Slater principle and it can be applied to both t-PA and NTs. It produces in t-PA the preferential direction in atom 
communication the only along chain axis (to be consequence of quasionedimensionality) , and it explains qualitatively 
the success of SSH-model. Since Slater principle can be equally applied to trans-polyacetylene (t-PA) chains and 
ZSCNTs, in the case of quasi-lD ZSCNTs it leads to applicability of SSH-model to given NTs, that is to NTs of small 
diameter despite the strong interaction between the chains. Experimental confirmation for given conclusion follows 
from ESR-studies in rather perfect CZSNTs, produced by high energy ion modifcation (HEIM) of diamond single 
crysatals 0], [Tlj ■ In the frames of SSH-model ZSCNT represents itself autonomous dynamical system with 

discrete circular symmetry consisting of finite number n £ N of carbon backbones of t-PA chains, which are placed 
periodically along transverse angle coordinate. Longitudinal axes {x{},i = l,n, of individual chains can be directed 
both along element of cylinder and along generatrix of any other smooth figure with axial symmetry. It is taken into 
account, that in the frames of SSH-model, the adjacent chains, which represent themselves a mirror of each other, 
will be indistinguishable, since the only one degree of freedom - dimerization coordinate u m of the m-th C-atom, 
m = I. N along chain molecular-symmetry axis x is substantial for determination of main physical properties in the 



2 



frame of SSH-model. 

The aim of given letter is to represent the experimental evidence for multichain QED-model, proposed in [l| and 
to confirm the theoretical prediction of Rabi wave phenomenon in 0]. The new phenomenon predicted is quantum 
coherent effect. It means, that along with requirement of quasionedimensionality the requirement of structural per- 
fectness arises. To observe quantum optical coherent effects on NTs, the ensemble of NTs has to be homogeneous. 
It means, that any dispersion in axis direction, chirality, length and especially in diameter both for single NT along 
its axis and between different NTs in ensemble has to be absent, axial symmetry has to be also retained, that is, 
there are additional requirements in comparison with, for example, t-PA technology. The CVD-technology of NTs 
production and many similar to its ones seem to be not satisfying to above-listed requirements at present. It means, 
that experimental results and their theoretical treatment will be different in both the cases, that really takes place 
(see for more details Sec. II). The situation seems to be analogous to some extent to the solid state physics of the 
same substance in single crystal and amorphous forms. The technology, based on HEIM, satisfy given requirements. 
Let us give short review concerning HEIM-method of production of incorporated carbon NTs. The first report on 
the discovery of new carbon phase - carbon nanotubes, incorporated in diamond matrix is related to 1990, and it was 
made during 1990 IBMM-Conference, Knoxwille, USA and it was repeated at E-MRS 1990 Fall Meeting, Strasbourg, 
France, [8J, that is, nothing was at that time known on the japan discovery of free NTs, related to 1991. It was found, 
that central place takes in ESR spectra of diamond single crystals, modified by high energy ion implantation, the 
strong line with very unusual radiospectroskopic properties. They was anisotropic, however the anisotropy was weak 
in comparison with the anisotropy of point centers in diamond. Especially interesting, that along with <?-value the 
linewidth AH pp was also found to be tensor quantitity, at that <?-tensor and A7? pp -tensor were characterized by the 
same symmetry group @. Moreover the only one eqivalent configuration was presented with principal direction along 
ion beam direction for both g and AH pp angular dependencies. Other lattice equivalent configurations were absent. 
The kind of symmetry group, corresponding to g-tensor and A_ff pp -tensor symmetry was strongly dependent on the 
choice of ion beam direction relatively the crystallographic axes of diamond lattice. It was found for the first time 
in radiospectroscopy, that g-tensor and AiJ pp -tensor symmetry are the mapping of the symmetry of lengthy objects 
with macrosizes along implantation direction [9] . Detailed studies of ESR spectral angular dependencies allowed to 
establish, that the structures formed by (lll)-HEIM in diamond have tracklike cylindrical shape, which are lengthy 
strongly in (11 Indirection, with the size varying from several mkm to several tens of mkm, depending on ion energy. 
In other words, the first studies in 1990 allowed to establish the formation of carbon nanotubes of cylindical symmetry 
shape incorporated in diamond matrix. Simultaneously the formation of quite different nanotubes was established. 
They represent themselves crimped cylinder with crimping corresponding to four-petal structure in the cross section. 
They are produced by the implantation direction, coinciding with (100) axis of diamond lattice The next step was 
the identification of the structure of cylindical symmetry shape (lll)-incorporated nanotubes with rolled up graphene 
sheet in zigzag shaped configuration. It was done by the study of radiospectroscopic properties of spin curriers in 
given tubes. It was found, that the system of paramagnetic centers (PC), which are responsible for appearance of 
strong absorption above described is non-Blochian system. A number of distinctive peculiarities have been observed 
for the first time in radiospectroscopy at all. The main ones among them are the following. 

Lit was found, that the shape of resonance lines cannot be described by Lorentzian, Gaussian, or by their con- 
volution. The shape is characterized by essentially more slow decrease of the absorption intensity on the wings in 
comparison with known shape functions. It was called super-Lorentzian. It was proved, that super-Lorentzian is gen- 
uine lineshape of resonance absorption (that is, it does not represent the superposition of Lorentzians with different 
lincwidths) [10| . 

2. It was established, that PC-system is nonsaturating with superlinear absorption kinetics (that is, the dependence 
of resonance signal amplitude on the amplitude of magnetic component of microwave field is superlinear) . Especially 
interesting, that superlinear absorption kinetics has been observed by both the registration of resonance signal in 
phase with high-frequency modulation and in quadrature with high-frequency modulation Q, [loj . 

3. It was found unusual for nonsaturating resonance systems dependence of resonance signal amplitude on the 
value of modulation frequency, which was characterized by substantial increase (instead decrease) of resonance signal 
amplitude with modulation frequency increasing Q, [Ioj |. 

4. It was established the effect of appearance of "phase angle", characterizing the absorption process. It consist 
in that, that despite nonsaturating absorption kinetics, the maximum of absorption is achieved not strictly in phase 
with modulation field, that is at modulation phase value, which is nonequal to zero. "Phase angle" is anisotropic in 
general case and maximum of its value was found to be equal 20° [Toj . 

The analysis of given peculiarities [TTJ] has lead to conclusion, that PC, which are responsible for strong 

non-Blochian absorption, are mobile and are characterized by two times of spin-lattice relaxation, by Ti, which is 
comparable with relaxation time of usual point PC in diamond single crystals (10~ 3 — 10~ 4 s) and by very short time 
r (upto 10 _13 s) of conversion of the energy of the spin system into mechanical kinetic energy of PC motion, which 
can be considered to be consequence of fundamental in magnetic resonance phenomenon of gyromagnetic coupling 
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between magnetic and mechanical moments. Therefore r-process is very fast process in comparison with any process, 
realized by means of pure magnetic interaction, the relaxation times of which are not shorter, than 10( - 10)s 0. 
The possibility of given conversion can be realized naturally, if the energy of the motion activation is very small, that 
is typical for mobile topological solitons like them identified in t-PA. So the conclusion, that the strong non-Blochian 
absorption is determined by mobile topological solitons was obtained by natural way. All subsecuent studies have 
confirmed given conclusion. It was found, that there is numerical coincidence of the characteristics of SSH topological 
7r-solitons in t-PA and in (111) incorporated NTs. So, the values of g-tensor components in Cu-implanted sample 
are g\ = 2.00255 (it is minimal g-value and it is g\ | principal direction of axial g-tensor, at that it coincides with 
ion beam direction), gi = g% = g± = 2.00273, the accuracy of relative g-value measurements is ±0.00002 |{|. g-value 
of paramagnetic 7r-solitons in trans-polyacetylene, equaled to 2.00263 [12], gets in the middle of given rather narrow 
interval of g-value variation of PC in ion produced NT. Although anisotropy of paramagnetic 7r-solitons in t-PA, 
which are also considered to be mobile paramagnetic centers, mapping the distribution of 7r-electron density along 
whole individual t-PA chain, is not resolved by ESR measurements directly (which in fact is the indication, that 
chemically produced t-PA is less perfect in comparison with NTs in diamond matrix), there are indirect evidences 
on axial symmetry of ESR absorption spectra in t-PA too, (l3j . fl4| . Consequently, the value 2.00263 is mean 
value and it coincides with accuracy 0.00002 with mean value of aforecited principal g-tensor values of PC in NTs. 
Given coincidence becomes to be understandable now, if to take into account the results of [l|, where is shown, that 
SSH-model is applicable to zigzag NTs, produced by HEIM of diamond single crystals. 

Thus we see, that it was possible to determine anisotropy of g- values with very high precision, that indicates on the 
very perfect and homogeneous NTs, produced with strict axial symmetry along ion beam direction and to establish 
the origin of paramagnetic centers in ion produced NT to be paramagnetic topological 7r-solitons (SSH-solitons), 
considered to be mobile PC, mapping the distribution of7r- electron density along whole individual NT [l(|, [llj . 
Thereby immediately from ESR studies was found, that structural element of new carbon phase with cylindrical 
symmetry produced in diamond by (111) high energy implantationis is t-PA chain backbone. At the same time t-PA 
chain backbone is basic element for both graphene, graphite, NTs formation and also for the carbyne formation. 
Carbynes were studied by ESR, they characterized by quite different ESR-spectra [HJ ■ The possibility of formation 
of cylindrical graphite rod is also not corresponding, since graphite is characterized by other ESR spectra. So it was 
concluded unambiguously, that the spectra observed by (111) high energy implantation are corresponding the only 
to graphene sheet, which is rolled up in that way, in order to t-PA chain backbones were directed along (lll)-axis of 
diamond lattice, coinciding with ion beam direction resulting in CZSNTs production. 

It is substantial, that free nanotubes were considered theoretically to be 2D-strutures and two dimensional lattice 
structure of a single wall carbon nanotube (SWNT) is specified uniquely by the chirality defined by two integers (n, m) 
[HI j 0- I n the Raman spectra of a SWNT arc two in-plane G point longitudinal and transverse optical phonon 
(LOandTO) modes [I?) and the out-of-plane radial breathing mode (RBM) [l8| are observed. The LO and TO 
phonon modes at the G point in the two-dimensional Brillouin zone are degenerate in graphite and graphene, however 
thy split in SWNT into two peaks, denoted by G + and G~peaks, respectively, [l9| because of the curvature effect. 
The agreement of experimental Raman studies of NTs with diameter > Inm with 2D SWNT-theory unambiguously 
indicates, that given NTs, produced by CVD are really 2D systems. At the same time the narrow perfect tubes with 
diameter < lnm cannot be considered strongly speaking to be 2Z?-systems, they are quasi-lD systems. Really, in [2(| 
is reported on the development of SWCNTs of 0.4 nm diameter - the smallest so far - inside the nanochannels of porous 
zeolite AlPOA — 5 single crystals and the authors are considered 0.4 nm diameter NTs to be ideal one-dimensional 
quantum hollow wires. The electronic and magnetic properties of these ultra-small NTs seem to be drastically different 
from those of large sized NTs. Really given tubes were shown to exhibit unusual novel phenomena like diamagnetism 
and superconductivity at low temperatures [21] in addition to several other optical properties and doping induced 
effects. It seems to be the consequence of onedimensionality of given NTs. The direct indication to given conclusion is 



recent ESR-studies in 22j , where electron spin resonance measurements on ultra-small single walled carbon nanotubes 



embedded in a SAPO 5 zeolite matrix with a main point of attention to potentially occurring CESR (ESR on the 
electrons in c-band) signals. Instead, only one paramagnetic signal is observed of symmetric shape at g = 2.0025 on 
the CNTs in zeolite cages, the ESR signal exhibits a predominantly Gaussian character at 8.9GHz(4.2K), with a 
AB, pp = 8.3G and g — 2.00251. For the free standing SWCNTs, the ESR signal takes a more Lorentzian shape at 
8.9GHz(4.2K ), now with AB pp = 5G and with g = 2.00254. So, from the above observed distinct changes in the 
ESR spectral properties of NTs with and without cages, it was inferred in [2^|, that the ESR signal indeed stems from 
the carbon NTs presenting inside the nanochannels of the zeolite. The observed g-value (2.00251) is inconsistent with 
the earlier measured values (g = 2.05 — 2.07) [23[ of CESR on non embedded single- and multiwalled CNTs of larger 
diameter. There was observed in [22j a dependence of peak-to-peak ESR signal width in confined SWCNT@SAPO 
5 at 70 if on microwave frequecy / suggesting a inhomogeneous (Gaussian) contribution to the line broadening. 
Least-square linear fit of ABpp is Bpp(G) — 7.7 + 0.018f (GHz). It is clear evidence, that the so called ultra-small 
carbon nanotubes, studied in [22j are not perfect. Nevetheless they can be really considered to be quasionedimensional 
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objects. It means, that like to t-PA, Peierls transition has to take place. Direct indication to given conclusion is g- 
value, which corresponds the only to deep centers in bandgap in any carbon materials [10|. On the one hand it is 
indicating of the appearance of bandgap in starting metal tubes to be consequence of Peierls metal-semiconductor 
transition. On the other hand g — 2.00254 is coincides in limits of accuracy of measurements) with g-value of g\ — 
2.00255 of axial g-tensor of SSH ir- -solitons in CZSNT, incorporated in diamond matrix by (lllWliamond HEIM and 
it near to the g- value of paramagnetic SSH 7r-solitons in trans-polyacetylene, equaled to 2.00263 [l2j, that indicates on 
the similarity of the structure of tubes incorporated in diamond and embedded in zeolite. At the same time angular 
dependences of g- value and AB pp were not observed in [22| . It seems to be additional argument to conclusion on 
insufficient homogeneity of NTs in cited work. Really they have different chiralities (authors report on the presence of 
three chiralities (5,0) (4,2) and (3,3)). On the other hand the characteristic for spin-glass temperature dependence of 
AB p p seems to be the indication, that the SSH tt -solitons are pinned and have random distribution along tube axis. 
For comparison, the temperature dependence of A.B pp in NTs, produced by HEIM is quite different and it indicates, 
that AH pp is not determined by unresolved g-anisotropy or unresolved hyperfine interaction [9(. It is interesting, that 
pinned solitons can be considered to be similar to chemical radicals, possessing by dangling bonds, which to some 
extent is in agreement with interpretation of PC-nature to be dangling bond defects in [22j, however they are not 
usual dangling bond point defects, but they are topological defects. Random distribution along tube axis means in 
turn, that individual tube withdefmed chirality is also not sufficiently perfect to observe optical quantum coherent 
effects. It means, that Slater principle with one preferential direction will be not applicable and optical spectra, in 
particular Raman spectra, will be determined the only by local properties of optical centers, that is, they can be 
similar to those ones in the tubes of larger diameters. It is substantial, that main Raman modes in the tubes of 
larger diameters seem to be similar in both perfect and in not sufficiently perfect NTs, since even in perfect NTs in 
given case two preferential direction are existing and onedimensional quantum coherent SSH-model |l| will be not 
applicable for description of optical properties in given tubes. 

II. RESULTS AND DISCUSSION 

Samples of type Ha natural diamond, implanted by high energy ions of copper (63 MeV, 5 x 10 14 cm~ 2 ) and boron 
(13,6 MeV) have been studied. Ion implantation was performed along (111) crystal direction. Raman scattering 
(RS) spectra were registered in backscattering geometry. Laser excitation wave length was 488 nm, rectangular slit 
350x350(/im) 2 was used, scan velocity was 100 cm -1 pro minute. 

The spectra observed are presented in Figures 1 to 3. All the spectra are 180° out of phase, that is dip positions 
correspond to maxima of the signals. The spectrum, presented in Figure 1, is characteristic of the only ion beam 
modified region (IBMR) of the sample, since the RS-line near 1332 cmT 1 , characteristic for diamond single crystals, 
is absent in the spectrum. The RS-lines with peak positions 656.8±0.2 cm -1 , 1215±1 cm -1 , 1779. 5±1 cm -1 and 
2022. 3±0. 5 cmT 1 correspond to optical centers in IBMR by laser excitation transversely to sample surface from 
implanted side. At the same time, by laser excitation of the same sample from opposite unimplanted side the 
qualitatively other picture is observed, see Figure 2. Firstly, the line, characteristic for diamond single crystals, is 
now presented in the spectrum (the whole amplitude of given line is not shown in Figure 2). It is interesting, that 
its frequency value (1328.7 cmT 1 ) is slightly shifted from usually observed value near 1332 cm -1 . It indicates on 
renormalization of optical phonon by Cu-implantation in the whole diamond matrix. 

Secondly, RS-lines with peak positions at 354.6, 641.8, 977.1 (±1 cm" 1 ), 1274.1 ±2 cm and more weak pronounced 
lines at 1569 ±3 cm -1 , 175 7±5 cmT x were observed. Moreover, all the lines listed are superimposed now with very 
broad (its linewidth is 1720±20 cmT 1 ) asymmetric line with peak position 1160±10cm _1 . The spectrum of boron 
implanted sample, Figure 3, was measured in the range (1000 - 2100) cm -1 . It is seen, that the spectra of boron 
and copper implanted samples are qualitatively similar in given spectral range. However numerical values - 1212. 3±1, 
1772. 5±1, 2011±0.5cm _1 - of the peak positions are slightly different, they are shifted to low frequency region. It is 
interesting, that there is regularity in peak position shift, which is increasing with frequency increase, and it is equal 
to 2.8, 7, 11.3 cmT 1 correspondingly. It is strong indication, that all three lines belong to the same optical system 
in both the samples studied. The line with peak position 1331.95 ±0.1cto _1 was also presenting. Its frequency value 
coincides with the value of optical phonon peak in conventional natural diamonds. The presence of given RS-linc 
by excitation from implanted side can be determined by two factors. Firstly, the effective thickness of IBMR is 
substantially less, it is ~ 1.4/im, in comparison with Cu-implanted sample (~ 7/im) @, consequently, laser excitation 
can reach an unimplanted region. It can also be suggested, that the modification is not entire in near surface region. 
Given results seem to be the first results in application of RS-spectroscopy to study of IBMR in diamonds. 

Comparing RS-spectra, represented in figures 1 to 3 with well known spectra of NTs produced by CVD-methods it is 
seen, that difference between the optical characteristics of the NTs, produced by HEIM of diamond single crystals and 
between the NTs, produced by other methods, including ultra-small NTs, is greatly despite the essential similarity of 
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FIG. 1: Spectral distribution of Raman scattering intensity in diamond single crystal, implanted by high energy copper ions, 
the excitation is from implanted side of the sample. 
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FIG. 2: Spectral distribution of Raman scattering intensity in diamond single crystal, implanted by high energy copper ions, 
the excitation is from unimplanted side of the sample. 

ESR-data with those ones for ultra-small NTs, embedded in zeolite matrix. On the other hand, comparing the figures 
1 and 2 between themselves, we can straight away conclude, that RS-lines observed cannot be attributed to usual 
atomic or molecular vibrations including the local vibration modes like to well known G-bands and radial breathing 
modes, observed in NTs, produced by CVD-methods, since the RS-spectra, corresponding to usual local vibrations 
cannot be dependent on the change of propagation direction of exciting light into opposite. Therefore 2D theory, 
developed for explanation of optical properties of NTs cannot be applicable for perfect NTs of small diameter. It 
is strong argument for application of the QED^model, proposed in [1[ for correct description of the objects studied. 
Since the theoretical model of small ID NTs [l[ is based on SSH-model of ID organic conjugated conductors let us 
consider SSH-model more carefully. It will be shown, that it needs in some corrections and development. 

We have repeated the calculation, presented in Q, Q. It has been found, that there is the second solution,which 
leads to independent branch of quasiparticles. Let us show it. We preserve all designation for quantities from [J]. We 
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FIG. 3: Spectral distribution of Raman scattering intensity in diamond single crystal, implanted by high energy boron ions, 
the excitation is from implanted side of the sample. 



have obtained the expression for Hamiltonian in reduced zone in the form 

= E 5>*( a * ~ + WfcA.pL - n v ks ) + 2NKu 2 e, (1) 



k 

where ctk and fiu are determined by 



Ek is determined by expression 



a k=\^A^k=\i^A- ( 2 ) 



in which 



Ek = J4 + ±l (3) 



£fc = 2to cos ka, Afe = 4cra sin ka (4) 



and n£ s , n£ s are occupation number operators in C- and F-bands correspondingly. The known SSH-solution corre- 
sponds to choice of lower signs in relations (j2J . The choice of upper signs gives independent solution. Therefore, for 
the energy of quasiparticles E^ u \ Efi u ^ in C- and V^-bands, corresponding to upper-sign-solution we obtain 



M 8EW(u) Al-ej v[u] 6EM(u) A| - c fc 



We see, that it differs from the energy of quasiparticles E^ l \ E^ in C- and F-bands, corresponding to lower-sign- 
solution (SSH-solution), which is 

k -~ML~- Ek ' E " ~ ~^KT " ~ ( ) 

E^ (u) and #M (u) in ([5]) and @ are eigenvalues of operator H(u), which correspond to upper and lower signs in ([2]). 
n ks ' n ks are eigenvalues of operators of particle numbers in C-band and V-band correspondingly. In [3| the coefficients 
Qffc and Pk were found from the conditions of minimum for the energy. However the only necessary conditions were 
used. At the same time the sufficient conditions for the minimum are substantial in given case, they change the role 
of both solutions. Really, since and /3k are coupled by the condition 

a 2 k +0l = l, (7) 
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we have conditional extremum and sufficient conditions for the minimum are obtained by standard study. The second 
differential of the energy to be the function of three variables a/-, @k and A, where A is coefficient in Lagrange 
function for conditional extremum, has to be positively defined quadratic form. From the condition of positiveness of 
three principal minors of quadratic form coefficients we obtain the following three sufficient conditions for the energy 
minimum 

a. The first condition 

{e k (l - §^)< ^|(nL - nl s ) < 0}, {e k (l - > ^\{n c ks - <J > 0} (8) 

&k &k &k &k 

for the SSH-solution and 



{e k (l + J < ^|(n c fes - nl s ) < 0}, {s k (l + f> |i|(n c fes - nt s ) > 0} (9) 

for the additional solution. It is seen, that the first condition is realizable for the quasiparticles of both the kinds, at 
that in equilibrium (n c ks — n v ks < 0) and in nonequilibrium (n ks — n ks > conditions. 

b. The second condition The second condition is the same for both the solutions and it is 

(4- - 2 ^) 2 - E 2 + - A A 2 k > (10) 

E k E k ' k A k K ' 

c. The third condition For the SSH-solution we have 

(3^+4-^)(nL-nL)>0. (11) 

It means, that SSH-solution is unapplicable for description of standard processes, passing near equilibrium state by 
any parameters. The quasiparticles, described by SSH-solution, can be created the only in strongly nonequilibrium 
state with inverse population of the levels in C- and F-bands. At the same time for the solution, which corresponds 
to upper signs in @, we obtain 

( 3 M_ 4 |f)(nL-nL)>0, (12) 

&k &k 

which is realizable both in near equilibrium and in strongly nonequilibrium states of the systems by corresponding 
choose of parameters. 

Let us consider the continuum limit for the ground state of the t-PA chain with quasiparticles of given branch. 
Taking into account, that in ground state n ks = 0, n ks = 1 we have 

E M( U ) = 7 A l7 £l dk + 2NKu\ (13) 

7T J ,/Af + el 



then, calculating the integral, we obtain 



E M {u) = *Mo {F{ * 1 Z 2) + i±4[£?(I, l-*)- F £, 1 - z 2 )]} + 2NKu\ (14) 

where 1 — z 2 ) is the complete elliptic integral of the first kind, -E(f , 1 — z 2 ) is the complete elliptic integral of 
the second kind, z 2 = ^j^- Approximation of (fT4| at z ^ 1 gives 

E M(u) =N{^- 6 -ln^^ + 2 ^ + ...} + 2NKu 2 . (15) 
7r 7r au to nto 

It is seen from (|15[) . that energy of quasiparticles, described by solution, which corresponds to upper signs in ^fy has 
the form of Coleman- Weinberg potential with two minima at the values of dimerization coordinate uq and — uq like 
to energy of quasiparticles, described by SSH-solution Q. 

It is remarkable, that SSH-model contains in implicit form along with the physical basis for the existence of solitons, 
polarons, breathers, formed in 7r-electronic subsystem, also the basis for the existence of similar quasiparticles in a- 
electronic subsystem. The origin is the same two-fold degeneration of ground state of the whole electronic system, 
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energy of which in ground state has the form of Coleman- Weinberg potential with two minima (|T5|) at the values of 
dimerization coordinate u and — u . Really the appearance of Uq ^ and — u Q ^ is concerned to the change of 
interatomic distance. It means, that simultaneously with n- subsystem, cr-subsystem is also dimerized. 

SSH-er-polarons have recently been found in carbynes [25| . and the formation of cr-polaron lattice (PL), which is 
antiferroelectrically ordered, was offered. Given chain state is characterized by the set of lines in IR-spectra, which 
were assigned with antiferroelectric spin wave resonance (AFESWR). Central mode is antiferroelectric resonance 
(AFR) mode, its value v°(C) in carbyne sample was equal to 477 cm -1 , the splitting parameter in IR detected 
AFESWR spectra was equal to 150 cm" 1 . Known values of v^(C) and AFESWR-splitting parameter in carbynes 
allow to estimate the range for expected values of analogous parameters in CZSNTs. AFESWR-splitting parameter 
is determined by exchange integrals in cr-electronic subsystem [2a ] . which seems to be the same in carbynes and 
in CZSNTs, since the role of quite different 7r-subsystems in carbynes and CZSNTs can be neglected to a first 
approximation. Consequently, the value of AFESWR-splitting parameter in CZSNTs (and in t-PA) has to be close to 
150 and to 300 cm" 1 by IR- and RS-AFESWR-detection correspondingly @. The frequencies i£(C) and v°{NT) 
of main AFESWR-mode in SSH-er-polaron lattice in carbyne and in CZSNTs depend on intracrystalline field |25j . 
that means, that their values will be different. However we can evaluate v°(NT), if to take into account the known 
relation for the vibration frequencies of similar centers, the fact of er-polaron and 7r-soliton lattice formation in 
carbynes, leading to change in effective masses [26j], presence of two 7r-subsystems in carbynes, difference of coherence 
lengths of cr-solitons and 7r-solitons (£>), band structure data for t-PA [29[ and carbynes [3(| ■ Relations for 
coherence lengths [3l[ are 

hv F hv F 

Cott = -r— ,4o CT = -t— , (16) 

where A 0(J , A 07r are a— and 7r-bandgap values at T = OK, vp is Fermi velocity. We obtain the frequency range 
for IR SSH-CT-polaron lines in t-PA and in CZSNTs v%{t - PA) € (386.7, 603) cm" 1 and v°(NT) G (402.5, 627.6) 
cm' 1 . Known IR-mode with the frequency near 540 cm" 1 , in t-PA [27j gets to interval (386.7, 603) cm" 1 and can 
represent itself AFR mode in cr-polaron lattice, that is, there is alternative interpretation of given IR-mode, ascribing 
earlier to Goldstone SSH-7r-soliton vibration mode [27j . From results of [HJ follows, that the same spectral interval 
with slightly different right-hand value, equaled to 673,7 cm" 1 , is the evaluation for the frequency of cr-polaron main 
AFESWR-mode in CZSNTs, which is RS-active. The calculation in [32| does not take into consideration the soliton 
and polaron formation. However cr-polaron formation does not violate the symmetry of task, that allows to conclude, 
that the same asymmetry value will be retained in IR and RS spectral distributions. Therefore, we come to conclusion, 
that intensive line 656.8±0.2 cm" 1 in RS-spectrum, presented in Fig.l, can be assigned with AFR mode of cr-polaron 
lattice, produced in NTs. The lines, 1215±1 cm" 1 , 1779. 5±1 cm" 1 and 2022. 3±0. 5 cm" 1 is revival part of Rabi 
wave packet in its frequency representation [l|. The arguments are the following. QED-Hamiltonian, proposed in [l[ 
for cr-polaron system of zigzag NT, is hypercomplex operator n-number 

[H] = [Ha] + [Hp] + [U*F] + \Uaa\- (17) 

The term [Ho] represents the operator of cr-polaron subsystem energy in the absence of interaction between cr-polarons 
themselves and with EM-field. The second term is the Hamiltonian of free EM-field, the third term describes the 
interaction of cr-polaron sybsystem with EM-field. The term [Hon] characterizes intrachain and interchain polaron- 
polaron interaction. Solution of the nonstationary hypercomplex Schrodinger equation for the state vector [l'F(i))] 
with given Hamiltonian in continuum limit can be represented in the form of sum of n solutions for n chains, where 
the solution for q-th chain & l q (x, t) is [l[ 

n-l 

$^(M) = ]>>^(M)[ei] p , (18) 
in which the matrix elements $> l qp (x,i) are determined by 

oo 

V77" flfil 

Q q (h, 0) exp expihxx 

(19) 

exp {i exp (dj (h) — g\Jl — Ikj (h))}dh, 
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where Q l q (h, 0),$j(/i), Kj(h) are determined by eigenvalues k Q G C, a = 0, n — 1 of n- numbers 
[ei] J is j-th power of the circulant matrix [ei] 



(h,0), 8{h) and x(h), 
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(20) 
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From equation ([3]) follows temporal dependence of the integral inversion for jth chain in NT and spectral dependence of 
the RS-signal amplitude, which explain the appearance of a number lines, which are additional to AFSWR lines in the 
spectra observed. The identification of the lines 1215±1 cm" 1 , 1779. 5±1 cm' 1 , 2022. 3±0. 5 cm" 1 (and corresponding 
lines in boron ion modified sample) with Fourier-image of revival part of Rabi packet is confirmed by the following. 
It is well known, that in the case of point absorbing centers Rabi frequency is linear function of the amplitude of 
oscillating EM- field. Given dependency takes also place for the center of Rabi packet, that follows from the analysis 
of Fourier transform of temporal dependence of the integral inversion. Further, it is evident, that amplitude of laser 
wave, penetrating in IBMR, is lesser by excitation from unimplanted side in consequence of some absorption in the 
unimplanted volume of diamond crystal. We see, that really, the values 1569 ±3 cm -1 , 1757±5 cm" 1 of two high 
frequency components are substantially less in given case, than 1779. 5±1 cm" 1 and 2022. 3±0. 5 cm" 1 , observed by 
excitation from implanted side, Fig. 2, Fig.l. Relative frequency changes are 1,151 and 1.134 ±0.003 for the pair 
[2022. 3±0. 5 cm -1 , 1757±5] and [1779. 5±1 cm" 1 , 1569 ±3] correspondingly and it is seen, that they are close to each 
other, at that larger frequency undergoes larger change in correspondence with theoretical analysis. The lines at 
354.6, 977.1 (±1 cm" 1 seem to be assigned with two AFESWR modes. The substantial decrease of relative intensity 
of 641.8 ±1 cm" 1 mode in comparison with 656.8±0.2 cm -1 mode testifies in favour of given assignment. It is seen, 
that AFESWR-splitting is rather large and it is comparable with the splitting between two polaron vibronic levels. 
It means, that linear AFESWR-theory [25|, which predicts a set of equidistant AFESWR-modes, arranged the left 
and the right of central mode, will be not applicable for given case. Really, AFESWR-modes are shifted on different 
distances 335.3, 287,2 cm -1 from main AFESWR-mode. At the same time average value of AFESWR-splitting is 
311.3cm -1 and it is close to the value in 300 cm -1 expected. Therefore, we obtain the direct proof of assigment of 
the lines 656.8±0.2 cm" 1 and 641.8 ±1 cm" 1 with cr-polaron lattice. 

To explain differences by the change of the direction of the excitation wave propagation, we have to take into 
consideration the following. Any SWR-excitation in ID-systems is strongly dependent on the directions of chain axis, 
vectors of intracrystalline external electric field, it is confirmed by ferromagnetic SWR study in carbynes, (24J. The 
axes in CZSNTs are not linear in the end of ion run and they are generatrixes of the figure of onion-like shape, 
that provides for necessary geometry for AFESWR-excitation of cr-polaron lattice. Moreover the appearance of very 
broad line seems to be indication on the excitation of the Frohlich movement of cr-polaron lattice itself. The presence 
of Frohlich sliding of cr-polaron lattice allows to explain qualitatively the appearance of "hysteresis" relatively the 
direction of excitation wave propagation from energy law conservation position. Really, moving cr-polaron lattice 
possesses by kinetic energy. It means, that it is required lesser energy value, to excite the local polaron vibration 
mode in correspondence with observation. 

Therefore CZSNTs represent themselves the example of the system, which strongly interact with EM-field. Ex- 
perimental detection of Rabi wave packets confirms on the one hand the theory, developed in [1| and it is also the 
first experimental confirmation of Rabi wave phenomenon, predicted in on the other hand. It means also, that 
semiclassical description of spectroscopic transitions in CZSNTs and in the systems like them cannot be appropriate, 
it substantially raises the practical concernment of QED-theory. The requirement of perfect onedimemensionality is 
substantial for applicability of SSH-model to CZSNTs, at the same time the Rabi wave phenomenon seems to be 
general and can be observed in 2D and 3D-systems. For instance, the first-order RS-peak in the vicinity of 1580 
cm -1 and so called second-order lines around 2480, 2700 (2D band) , 3250 cm" 1 by \ exc — 488 nm were observed 
in graphene [33[ , at that the intensity of 2D band exceeds substantially the intensty of the first-order line. Moreover 
theoretical dependence of 2D band position on excitation energy, presented in [33[ (Fig. 5), is far from experimental 
dependence. At the same time the suggestion, that the second-order lines are spectral mapping of revival part of Rabi 
wave packet is agreeing with data in Fig. 5 very well, since, if to choose instead energy the electrical component of 
oscillating excitation-field to be x-coordinate, then we obtain linear experimental dependence in correspondence with 
expected for the center of Rabi packets. Similar explanation can be proposed for transitions in heavily boron-doped 
diamond in the region of 1400-2800 cm -1 (34J. Naturally Rabi waves can also be observed in CVD-NTs. 

Thus, QED-model for multichain coupled qubit quasi-lD system, proposed in [l[, is confirmed by Raman scattering 
studies of carbon zigzag-shaped nanotubes, produced by HEIM of natural diamond single crystals. New quantum 



10 



optics phenomenon - Rabi waves, predicted in [2|, has been experimentally identified for the first time. 
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